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Abstract 

Demands from several areas of science and technology have lead to a worldwide search for accurate 
optical clocks with an uncertainty of 1 part in 10^*, which is 10'' times more accurate than the 
present day cesium atomic clocks based on microwave frequency regime. In this article we discuss 
the electric quadrupole and the hyperfine shifts in the 5s^Si/2 — > 4d ^1)5/2 clock transition in Sr""", 
one of the most promising candidates for next generation optical clocks. We have applied relativistic 
coupled cluster theory for determining the electric quadrupole moment of the 4d ^1)5/2 state of ^^Sr"*" 
and the magnetic dipole {A) and electric quadrupole {B) hyperfine constants for the 5s^5'i/2 and 
states which are important in the study of frequency standards with Sr"^. The effects of 
electron correlation which are very crucial for the accurate determination of these quantities have 
been discussed. 

PACS number(s). : 31.15.Ar, 31.15.Dv, 32.30.Jc, 31.25.Jf, 32.10.Fn 

1 Introduction 

The frequencies at which atoms emit or absorb electro-magnetic radiation during a transition can be 
used for defining the basic unit of time [11 O [3]. The transitions that are extremely stable, accurately 
measurable and reproducible can serve as excellent frequency standards [U [2] . The current frequency 
standard is based on the ground state hyperfine transition in ^^'^Cs which is in the microwave regime 
and has an uncertainty of one part in 10^^ [1]. However, there is a search for even more accurate clocks 
in the optical regime. The uncertainty of these clocks is expected to be about 1 part in 10^'' or 10^^ 
[5]. Some of the prominent candidates that belong to this category are ^^Sr"'" [6l[7], ^^^Hg+ [8], ^^^Yb+ 
[9], ^^Ca+ [To], ^'^^Ba+ etc. Indeed detailed studies on these ions will have to be carried out in order to 
determine their suitability for optical frequency standards. In a recent article [H] Gill and MargoHs have 
discussed the merits of choosing ®^Sr+ as a candidate for an optical clock. Till recently, the most accurate 
measurement of an optical frequency was for the clock transition in ®^Sr+ which has an uncertainty of 
3.4 parts in 10^^ [12j. However, recently, Oskay et al. [13] have measured the optical frequency of ^^^Hg"'' 
to an accuracy of 1.5 parts in 10^^ and further improvements are expected [14j . 

In this article we concentrate on strontium ion (Sr+) which is considered to be one of the leading 
candidates for an ultra high precision optical clock [TT]. The clock transition in this case is ?>s'^Si/2 — > 
4^^1)5/2 and is observed by using the quantum jump technique in single trapped strontium ion. When 
an atom interacts with an external field, the standard frequency may be shifted from the resonant 
frequency. The quality of the frequency standard depends upon the accurate and precise measurement 
of this shift. To minimize or maintain any shift of the clock frequency, the interaction of the atom with 
it's surroundings must be controlled. Hence, it is important to have a good knowledge of these shifts 
so as to minimize them while setting up the frequency standard. Some of these shifts are the linear 
Zeeman shift, quadratic Zeeman shift, second-order Stark shift, hyperfine shift and electric quadrupole 
shift. The largest source of uncertainty in frequency shift arises from the electric quadrupole shift of the 
clock transition because of the interaction of atomic electric quadrupole moment with the gradient of 
external electric field. In this article we have applied relativistic coupled-cluster (RCC) theory, one of 
the most accurate atomic many-body theories to calculate the electric quadrupole moment (EQM) and 
the hyperfine constants for the energy levels involved in the clock transition. 
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Figure 1: Clock transition in ®^Sr+. Excitation energies of the 4d ^I?3/2and id ^Z35/2levels are given in 
cm~^. 

2 Relativistic coupled-cluster theory 

The relativistic and dynamical electron correlation effects can be incorporated in many-electron systems 
through a variety of many-body methods [iSl [TBI [TTl [18] . The relativistic coupled cluster (RCC) method 
has emerged as one of the most powerful and effective tools for high precision description of electron 
correlations in many-electron systems \T7\ 118). The RCC is an all-order non-perturbative scheme, and 
therefore, the higher order electron correlation effects can be incorporated more efficiently than using the 
order- by-order diagrammatic many-body perturbation theory (MBPT). RCC is equivalent to all order 
relativistic MBPT (RMBPT). The RCC results can therefore be improved by adding the important 
omitted diagrams with the aid of low order RMBPT. We have applied RCC theory to calculate atomic 
properties for several systems and more details can be obtained from Ref. [19] . 

Here we present a brief outline of RCC theory. We begin with A^-electron Dirac-Coulomb Hamiltonian 
(H) which is expressed as 

N N n 



H = ^[cdi-p, + f3mc^ ^VN{,ri)\+^ — , (1) 



1=1 

with the Fermi vacuum described by the four component Dirac-Fock (DF) state |<I>). The normal ordered 
form of the above Hamiltonian is given by 

H^^H- imi'^) = E(*I/IJ) + i ^ {ala]aia,} , (2) 

where 

mki) = (tj\—\ki}-{ij\—\ik}. (3) 

Following Lindgren's formulation of open-shell CC [20], we express the valence universal wave operator 
as 

n = {exp(a)}, (4) 

and a being the excitation operator and curly brackets denote the normal ordering. The wave operator 
fl acting on the DF reference state gives the exact correlated state. The operator a has two parts, one 
corresponding to the core and the other to the valence sector denoted by T and S respectively. In the 
singles and double (SD) excitation approximation the excitation operator for the core sector is given by 

T = T,+T2 = J2 {«P«a} € + H'^l'^baa} til , (5) 

ap abpq 

tP and t^l being the amplitudes corresponding to single and double excitations respectively. For a single 
valence system Hke Sr+, the excitation operator for the valence sector turns out to be exp(S') = {1 + 5*} 
and 
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Figure 2: Some many-body diagrams representing the electric quadrupole/hyperfine interaction. Holes 
(occupied orbitals, labeled by a) and particles are denoted by the lines directed downward and upward 
respectively. The double line represents the interaction vertex. The valence (labeled by v) and virtual 
orbitals (labeled by p,q,r..) are depicted by double arrow and single arrow respectively whereas the 
orbitals denoted by © can either be valance or virtual. 



S = Si + 82^ 22 + zl {aWi^bak} sll , (6) 

k^p bpq 

where and s^l denotes the single and double excitation amplitudes for the valence sectors respectively. 
In Eqs. ^ and ^ we denote the core (virtual ) orbitals by a,b,c... {p,q,r...) respectively and k corre- 
sponds to the valence orbital. 

The corresponding correlated closed shell state is then 

Ivf) = exp(r) 1$) . (7) 

The exact open shell reference state is achieved by using the techniques of electron attachment. In order 
to add an electron to the kth virtual orbital of the N electron DF reference state we define 

i*r'>^4i'i') (8) 

with the particle creation operator a\. Then by using the excitation operators for both the core and 
valence electron the (A^ 1) electron exact state is defined as [20] : 

K+i)=exp(r){l + 5}|<i>f+i). (9) 

We write the expectation value of any operator O in a normahzed form with respect to the exact state 
1*^+1) as 

(^w+i|0|^w+i) _ {1 + 5t} exp(Tt)Oexp(T) {1 + S} 

\ ^ ~ (^.Jv+i| ^jv+i^ - {1 + S-t} exp(rt) exp(r) {1 -f 5} ' 

For computational simplicity we store only the one-body matrix element of O = exp(T^)0 exp(T). O 
may be expressed in terms of uncontracted single-particle lines j2T]. The fully contracted part of O will 
not contribute as it cannot be linked with the remaining part of the numerator of the above equation. 
First few terms of Eq. l|10p can be identified as O , OSi and OS2 which we identify as dressed Dirac- 
Fock (DDF), dressed pair-correlation (DPC) and dressed core-polarization (DCP) respectively. The other 
terms are being identified as , S'jOS'i, S'jOS'2, SlOSi and S\0S2 and are classified as higher order effects. 
We use the term 'dressed' to describe the modification of the operator O due to core-core correlation 
effects. In Fig. [2]we replace the operator O by the dressed operator O (the dressed quadrupole/hyperfine 
interaction operator) which includes the core excitation effects and the respective figures are termed as 
Ha) OSi,Wih+c) OS2 and Eld) SlOSi. Here Ha) and Hb+c) represent the DPC and DCP effects 
respectively. [Hb) is known as the direct and^c) is the exchange DPC diagram. Figure Hd) refers to 
one of the higher order pair correlation effect which belongs to the set, termed as 'others'. 

3 Electric quadrupole shift 

The largest source of systematic frequency shift for the clock transition in Sr+ arises from the electric 
quadrupole shift of the 4:d'^D^/2 state caused by its electric quadrupole moment of that state and the 
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interaction of the external electric field gradient present at the position of the ion. The electric quadrupole 
moment in the state Ad ^1)5/2 was measured experimentally by Barwood et al. at NPL [22j. Since 
the ground state 5s ^S'1/2 does not possess any electric quadrupole moment, the contribution to the 
quadrupole shift for the clock frequency comes only from the 4d state. 

The interaction of the atomic quadrupole moment with the external electric-field gradient is analogous 
to the interaction of a nuclear quadrupole moment with the electric fields generated by the atomic 
electrons inside the nucleus. In the presence of the electric field, this gives rise to an energy shift by 
coupling with the gradient of the electric field. Thus the treatment of electric quadrupole moment is 
analogous to its nuclear counterpart. The quadrupole moment of an atomic state |^'(7, J, M)} is 
defined as the diagonal matrix element of the quadrupole operator with the maximum value Mj, given 
by 



© = (*(7jj)|e,,|*(7jj)) . (11) 

Here 7 is an additional quantum number which distinguishes the initial and final states. The electric 
quadrupole operator in terms of the electronic coordinates is given by 

3 

where the sum is over all the electrons and z is the coordinate of the jth electron. To calculate the 
quantity we express the quadrupole operator in its single particle form as 

el^^^E^i^ (13) 

m 

More details about evaluation of electric quadrupole moment using RCC theory is described in our recent 
paper [23]. The electric quadrupole shift is evaluated using the relation 



[{2F + 3){2F + 2){2F + 1)2F(2F - 1)]^^^ 

and 

0{a, (3) = [(3 cos^ /3 - 1) - e{cos^ a - sin^ a)] . (15) 

Here 7 specifies the electronic configuration of the atoms and F and Mp are the total atomic angular 
momentum (nuclear + electronic) and its projection; a and (3 are the two of the three Euler angles that 
take the principal-axis frame of the electric field gradient to the quantization axis and e is an asymmetry 
parameter of the electric potential function [2i] . 



4 Hyperfine shift 

The frequency standard is based on the *^Sr isotope. In addition to ^^Sr"*", the odd isotope ^^Sr+ has also 
been proposed as a possible candidate for an optical frequency standard [25|- An experiment has been 
performed in NPL to measure the hyperfine structure of the Ad state in ®'^Sr+[25]. Theoretical 

determination of hyperfine constants is one of most stringent tests of accuracy of the atomic wave 
functions near the nucleus. Also accurate predictions of hyperfine coupling constants require a precise 
incorporation of relativistic and correlation effects. 

Unlike the even isotope (^^Sr"*") of strontium ion, ^^Sr"*" has a non zero nuclear spin (/ = |) and the 
mp = levels for both the ^5*1/2 and ^£>5/2 states are independent of the first order Zeeman shift. Here, 
in table [2] we present the results of our calculation of the magnetic dipole {A) hyperfine constant for the 
5s^S'i/2 and Ad ^1)5/2 states and the electric quadrupole hyperfine constant {B) for the Ad state 
of ®^Sr+ and compare with the measured values. More details of our calculation can be found in Ref 

The hyperfine interaction is given by 

77;,/, = ^mW-TW, (16) 

k 
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where M^^) and T^*^) are spherical tensors of rank k, which corresponds to nuclear and electronic parts 
of the interaction respectively. The lowest fc = order represents the interaction of the electron with the 
spherical part of the nuclear charge distribution. 

In the first order perturbation theory, the energy corresponding to the hyperfine interaction of the 
fine structure state \JMj) are the expectation values of H/ifg such that 



W{J) = {IJFMf\ ^Aff'^) ■ T^'^) \IJFMf) 



k 

^(-1)^+'^+^^ '\ I [(/||AfW||/)(J||TW||J) 



(17) 



Here I and J are the total angular angular momentum for the nucleus and the electron state, respectively, 
and F = I + J with the projection Mp- 

4.1 Magnetic dipole hyperfine constant 

Hyperfine effects arise due to the interaction between the various moments of the nucleus and the 
electrons of an atom. Nuclear spin gives rise to a nuclear magnetic dipole moment which interacts with 
the electrons and thus gives rise to magnetic dipole hyperfine interaction defined by the magnetic dipole 
hyperfine constant A. For an eigen state \IJ) of the Dirac-Coulomb Hamiltonian, A is defined as 



V / / v/J(J+l)(2J+l)' ^ ' 



where /x/ is the nuclear dipole moment defined in units of Bohr magneton /ijv. The magnetic dipole 
hyperfine operator Tq^^ which is a tensor of rank 1 can be expressed in terms of single particle rank 1 
tensor operators and is given by the first order term of Eq. lfT7|) 

^i^^ = E4^' = E— (19) 

9 3 

Here 7? is the Dirac matrix and Y^^ is the vector spherical harmonics. The index j refers to the j-th 
electron of the atom with Vj its radial distance and e is the magnitude of the electronic charge. 

4.2 Electric quadrupole hyperfine constant 

The second order term in the hyperfine interaction is the electric quadrupole part. The electric quadrupole 
hyperfine constant is defined by putting fc = 2 in Eq. (fTT]) . The nuclear quadrupole moment is defined 
as 

-Ef^^E— ('^o, (20) 



Here, Cg'^^ — \J (2k+i) ^kq^ with Y^q being the spherical harmonic. Hence the electric quadrupole hyperfine 
constant B in terms of the nuclear quadrupole moment Qn is 

2J(2J-1) "1^^^ 



The corresponding shift in the energy levels are known an hyperfine shift which is expressed as 

Wuy, - Wmi + We2 - a- + 2/(2/-l)2J(2J-l) ' ^^^^ 

where K = F{F + 1) - /(/ + 1) - J(J + 1). 
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Table 1: Electric quadrupole moment for the id "^D^ji state of *^Sr+ in units of ea§. PW corresponds 
to our present CCSD(T) calculation and MCDF for Multi-configuration Dirac-Fock. 

FW MCDF [2f| Experiment [22j 
4^5/2 2^94 3^02 2.6±0.3 



5 Results and discussions 

The occupied and the virtual orbitals used in the calculation are obtained by solving the Dirac-Fock 
(DF) equation for Sr++ for a finite Fermi nuclear distribution. These orbitals are linear combinations 
of Gaussian type functions on a grid [28j. The open shell coupled cluster (OSCC) method is used to 
construct different single valence states whose reference states correspond to adding a particle to the 
closed shell reference state. We use the singles-doubles and partial triples approximation, abbreviated 
as CCSD(T) and excitations from all the core orbitals have been considered. We have estimated the 
error incurred in our present work, by taking the difference between our RCC calculations with singles, 
doubles as well as the most important triple excitations (CCSD(T)) and only single and double excitations 
(CCSD). 

5.1 Electric quadrupole moment 

We present our results of the electric quadrupole moment for the \d state of ®*Sr+ in table [TJ 

The value of in the Ad state measured experimentally is (2.6 ± 0.3)eaQ [22]. Our calculated 

value for the Ad stretched state is (2.94 ± 0.07)eaQ, where e is the electronic charge and ao is the 

Bohr radius. We analysed our results and have found that the DDF contribution is the largest. The 
leading contribution to electron correlation comes from the DPC effects and the DCP effects are an order 
of magnitude smaller. This can be understood from the DPC diagram (Figl^a)) which has a valence 
electron in the 4^5/2 state. Hence the dominant contribution to the electric quadrupole moment arises 
from the overlap between virtual d^/2 orbitals and the valence, owing to the fact that Si is an operator 
of rank and the electric quadrupole matrix elements for the valence 4^5 /2 and the diffuse virtual d^ /2 
orbitals are substantial. On the other hand, in the DCP diagram (Figl2]Jb+c)), the matrix element of 
the same operator could also involve the less diffuse s or p orbitals. Hence, for a property like the electric 
quadrupole moment, whose magnitude depends on the square of the radial distance from the nucleus, 
this trend seems reasonable, whereas for properties like hyperfine interaction which is sensitive to the 
near nuclear region, the trend is just the opposite for the d states |29]. As expected, the contribution 
of the DHOPC effect i.e., S\OSi (Figl2]Jd)) is relatively important as it involves an electric quadrupole 
matrix element between the valence 4^5/2 and a virtual ^5/2 orbital. 

5.2 Hyperfine constants 

The magnetic dipole {A) and electric quadrupole [B) hyperfine constants for the of *^Sr+ are given 
in table [2] along with the calculated and experimental results. The gi = fiN {^) value used for the 
calculation is from Ref. [30]. PW corresponds to the 'present work' using CCSD(T). CC stands for 
coupled-cluster calculation by Martensson-Pendrill [3T| and Nayak and Chaudhuri f32], DF-AO for Dirac- 
Fock with all order core polarization effect by Yu et al. and the column 'others' refers to the calculation 
by Yuan et al. using relativistic linked cluster many-body perturbation theory (RLCMBPT) [S^ and 
by one of the authors using relativistic effective valence shell Hamiltonian method |34] . Information on 
A for the Ad ^D^/2 state is very important in connection with optical frequency standards |12[ [25]. The 
measured value of this quantity is 2.1743 ± 0.0014 MHz [25], whereas the previously calculated values 
vary from 1.07 MHz [ST] and 2.507 MHz [35]. Our calculated value of A is 2.16 ± 0.02 MHz; this is the 
most accurate theoretical determination of A for the Ad ^D^/2 state to date. For the 5s ^S'i/2State our 
calculated value of A is 997.26 ± 0.03 MHz. To analyze the result we focus on the various many-body 
effects contributing to the calculation of A. The most important many-body diagrams are presented 
in Fig. [2I We have noticed that for 5s ^<S'i/2 state the dominant contribution is at the DDF level 
~ 72%. However, for the Ad ^D^/2 state the DCP effect is larger and its sign is opposite that of the 
DDF contribution. In their calculations Martensson [31] and Yu et al. [35] have found similar trends. 
We have seen that the higher order correlation effects contribute significantly in determining A for this 
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Table 2: Magnetic dipole (A) hyperfine constant for the 55^5*1/2 and Ad states and the electric 

quadrupole hyperfine constant {B) for 4(i ^£>5/2 state of ®^Sr+ in MHz. 



4^5/2 A 



B 



PW 



CC DF-AO |35j Others 



Experiment 



5si/2 A -997.26 



2.16 
47.8 



-1000 [31] 
-999.89 132] 



1.07 [31] 
1.87 [32] 

54.4 [m, 

51.12 [32] 



-1003.18 



2.51 



-987 [33] 
-1005.74 134J 



-iooo.5±i.o [sg 

-990 [37] 
-993.5 [38] 
-1000.473 673 [39] 

2.1743±0.0014 [25j 

49.11±0.06 [361 



state - collectively they are 60% of the total value but opposite in sign. In the earlier calculations the 
determination of the higher order effects was not as accurate as ours. 

The calculated value for the electric quadrupole hyperfine constant {B) for the Ad ^1)5/2 state is 
47.8 ± 0.2 MHz which deviates ~ 2.7% from the central experimental value. The earlier determination 
of B was off by ~ 11% from the experiment. Since the other state involved in the clock transition is 
spherically symmetric the electric quadrupole hyperfine constant is zero and there will be no hyperfine 
shift due to B for the 5s ^S'i/2State. 

From the figure presented in Fig. [2] it is clear that the DPC effect involves the hyperfine interaction 
of a valence electron and the residual Coulomb interaction, i.e for 5s ^S'1/2 state the hyperfine matrix 
element becomes ([4p^]5si/2 |/i/i/s| [4p^]9), where q can be any virtual orbital and h^fs is the single- 
particle hyperfine operator. Since only Si/2 and pi/2 electrons have a sizable density in the nuclear 
region, the DPC effect is dominant for 5si/2 state but not for the 4^5/2 state. On the other hand, the 
DCP effects represent the hyperfine interaction of a polarized core electron with any virtual electron 
(see Fig. [21d,c). For 4^5/2 state it is clear that the hyperfine matrix element for the DPC effect is much 
smaller than the DCP effect, which plays the most important role in determining the value of A. The 
core polarization contribution is so large for this state that it even dominates over the DDF contribution. 

6 Conclusion 

In conclusion, we have performed an ah initio calculation of the electric quadrupole moment for the 
Ad state of ®*Sr+ to an accuracy of less than 2.5% using the RCC theory. This is the first appHcation 

of RCC theory to determine the electric quadrupole moment (EQM) of any atom and is currently the 
most accurate determination of EQM for the Ad ^-05/2 state in Sr"*". We have also calculated the magnetic 
dipole {A) hyperfine constant for the 5s ^S'1/2 and Ad ^£>5/2 states and electric quadrupole hyperfine 
constant {B) for the Ad ^B)^i2 state of *^Sr+. Evaluation of correlation effects to all orders as well as the 
inclusion of the dominant triple excitations in our calculation was crucial in achieving this accuracy. The 
magnitude of electric quadrupole moment depends on the square of the radial distance from the nucleus, 
whereas properties like hyperfine interaction are sensitive to the near nuclear region. The accurate 
determination of quantities like electric quadrupole moment and hyperfine constants establish the fact 
that RCC is very powerful and efficient method for determining atomic properties near the nuclear region 
as well as at large distances from the nucleus. Our result will lead to a better quantitative understanding 
of the electric quadrupole shift of the resonance frequency of the clock transition in Sr+ . 

Acknowledgments : Financial support from the BRNS, DAE for project no. 2002/37/12/BRNS 
is gratefully acknowledged. The computations are carried out in our group's Xeon and 
the Opteron computing cluster at IIA. 
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